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Abstract 

We put disks on a sphere between two parallels of this sphere. The disks have the curvature of 
the sphere and interact via a simplified Hertz law. We analyze the behavior of the total kinetic 
energy of the whole assembly of disks with Riemanian geometry and contrarily to flat spaces, this 
kinetic energy increases. This is due to a non radial component of the resulting velocity vector at 
each interaction between two disks. The study of granular matter in twodimensional curved spaces 
is therefore not trivial. 

PACS numbers: 45.70.-n;02.40.Ky;45.20.dh 
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Granular matter has been a wide subject of interest since the 80 's |l|,y,|3l. But up to 
now no study of granular matter and in particular soft particles has been made in a space 
with curvature. The rare studies which deal with curved spaces are the articles of Caillol ^ 
for plasmas on spheres and hyperspheres. Although, these studies are very different from 
granular matter analysis as electrons have long range interactions. Therefore the fact that 
the space are spheres or hyperspheres is linked to boundary conditions j^. 

For granular matter, the interactions are short ranged and even if one tried to compare 
periodic boundary conditions for a plane with a sphere, the curvature interferes and the 
comparison can not be direct (essentially because of the metric). 

Here we used a Molecular Dynamics program for soft disks on a sphere. A simplified 
Hertz law has been used in order to match the metric of differential geometry of the sphere. 
In order to check that our results were not an artifact of the computational program we used 
two algorithm for the integration: the predictor corrector and the velocity Verlet algorithms. 
Both gave the same results. 

We put disks with equal radius on spheres with different radii. The spheres match 
exactly the curvature of the sphere, i.e. the space where they are enclosed is the sphere. 
The number of spheres may vary. The surface of the spheres where they can move is limited 
by two parallels of the spheres at equal distance from the equator, and that in order two 
avoid the poles which are critical points for the calculations (see fig.l). Gravity acts from 
north pole on the top of the sphere to south pole on the bottom of the sphere. The disks 
are soft and interact via a simplified Hertz law in a curved space. 

In order to take account of the curvature, we used the following metric tensor of Riemanian 
geometry: 

Tn = 1 (1) 
T22 = (2) 
n^ = r^sm^e (3) 

The other components of the tensor are equal to zero. For all physical parameters included 
in the following equations and using lengths dimensions, we used this tensor. 

The simplified Hertz law is written also as follows, always in the spirit of the metric 
tensor: 

Fij = [Kn{d-rij) -me7„(vij.n)]n + emm[7^|vij.s|,/i||F„|]s (4) 



with Fij is the relative force between disks i and j, Kn is a spring constant, d is the diameter 
of the disks, rij is the center to center distance between disks i and j, nie is the mass of 
one disk, 7„ (resp. 7^) is the dissipation coefficient in the normal ( resp. shear)direction, 
n (resp. s) is the normal (resp. shear) vector (with respect to two disks), and Fn is the 
normal force. 

e = (5) 
|vij.s| 

e is the sign of the relative velocity in the shear direction s. 

The two parallels arc located at 6 = 71/ 2 and 6 = 37r/2. The disks interact with the 
parallels with a damping coefficient equal to 0.75. The data are have no dimensions taking 
account of their mass m, the gravity g and their diameter d. So, — 2000'mg/d, 7^ = 
25{d/gy/^,js = '^5{d/gy/^ and n = 0.45. 

The normal vector is calculated as follows: 



di rdO ^ rsini^vx. , . , 

n = d£ = — = —d + ^^(/) 6) 

' \di\ \di\ 



and the shear vector is calculated as follows: 

r sin Odd) ^ rd6 , 
s = ^—-9 —d) (7) 

\di\ \di\ 

therefore n.s = in spherical coordinates r, 9 and 0. 

At the beginning of computation, i.e. at Molecular Dynamics Step (MDS) equal to 1, the 
disks are put at random between the two parallels with a random distribution of velocity 
vectors in all directions. Gravity is acting since the beginning of computation. 

When the disks are put at random within the two parallels and do not see each other, 
we checked that they fall to the bottom parallel until reaching it. We also checked that two 
disks which have a trajectory exactly on the equator with the same velocity vector norm 
but with inverse direction, the acceleration of the two disks did fall to zero as well as their 
kinetic energy. With the same program, but with a radius of the sphere infinite (a plane), 
we checked that the kinetic energy deduced from the acceleration of the disks was tending 
to zero as is usually find in a plane. 

The results are shown in figures 2 and 3. Figure 2 corresponds to different radii of the 
sphere but with a constant number of disks. As one can see, for increasing radius of the 
sphere, the maximum of the kinetic energy is lower than for a smaller radius of the sphere. 
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Figure 3 corresponds to different number of disks but with an equal radius of the sphere. 
Here again, for an increasing number of disks the resulting maximum kinetic energy is larger. 

Let us have an explanation of this phenomenon of increasing kinetic energy (not tending 
to zero though the existence of dissipation in the Hertz law). For that, we plotted figure 
4 which represents two disks in interaction on the surface of the sphere but with different 
norms of the velocity vector after interaction. If one calculates the resulting vector of the 
two interacting velocity vectors after interaction, one can see that this resulting vector has a 
non radial component with respect to the center of the sphere. This means that at after each 
interaction, the non radial component obtained will give place to a tangential component of 
the resulting velocity vector of each disk. This fact will give an additional term to the total 
kinetic energy of the whole assembly of disks. Therefore the kinetic energy will increase 
infinitely. In our case, as we use the Hertz law with dissipation, the kinetic energy will 
stabilize at a value where the number of interactions and the value of the additive kinetic 
energy for each time step will be counterbalanced by the dissipation. 

In this letter, we showed that unlike plane spaces, in a curved space which is a sphere, 
the total kinetic energy of a whole assembly of disks which are on the sphere and have the 
curvature of the sphere increases until reaching a stable value in time. This value is different 
from zero and is higher for larger numbers of disks and smaller sphere radii. The explanation 
of this phenomenon is the non radial component which appears at each interaction following 
the Hertz law, between two disks, and which leads to an additive non zero kinetic energy to 
the total kinetic energy. Studying twodimensional disks in a space with non zero curvature 
is therefore not trivial and very different from the study of granular matter in flat spaces. 

Acknowledgement We would like to acknowledge financial support for the PCP Statique 
et dynamique des milieux granulaires 



4 



[1] L. P. KadanofF, Review of Modern Physics, 71, 435, (1999); and references within 

[2] S. Nagel, Review of Modern Physics, 64, 321, (1992); and references within 

[3] H. M. Jeager, S. R. Nagel and R. P. Behringer, Physics Today, 49, 32, (1996); Review of 

Moderm Physics, 68, 1259, (1996); and references within 
[4] J. M. Caillol, J. Chem. Phys. Ill 6528 (1999) 
[5] J. M. Caillol and D.Gilles, J. of Stat. Phys. 100 905 (2000) 



5 




FIG. 1: Schematic representation of the disks on the sphere hmited by two parallels. The disks 

are represented as small spheres because of software displaying 
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FIG. 2: Total kinetic energy as a function of time for the different sphere radii (in unit of disk 
diameter) but a constant number of disks 
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FIG. 3: Total kinetic energy as a function of time for a constant sphere radius (in unit of disk 
diameter) but different number of disks 
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FIG. 4: Schematic representation of the non radial component of the relative velocity vector after 
an interaction between two disks on the sphere 
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